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§ I ABSTRACT 

Of the known pulsar wind nebulae, 8 are good candidates for being in the early stage of evolution where the wind 

\ nebula is interacting with the freely expanding supernova ejecta. Several of these have been identified with historical 

■ supernovae. Although the identification of SN 1181 with 3C 58 has been thought to be relatively secure, the large size 

O of the nebula, the amount of swept up mass, and the internal energy indicate a larger age. For Gl 1.2-0.3, the nebular 

^ | size and internal energy are consistent with the identification with the possible supernova of 386. Although the Crab 

y— ( ■ Nebula appears to have approximate energy equipartition between particles and the magnetic field, the nebulae 3C 

O , 58 and MSH 15-52 appear to be particle dominated. The low magnetic field is consistent with models in which the 

^ | nebulae are created by a shocked pulsar wind. 

^: 

O INTRODUCTION 

Pulsars are expected to be born inside the s upernova explosions of massive stars, which provide the surroundings 



X 



for the initial evolution of their wind nebulae jRevnolds and Chevalier. 1984 ). The PWNe (pulsar wind nebulae) 



initially expand in the freely expanding ejecta of the supernova. Eventually, the reverse shock wave from the supernova 
interaction with the surrounding medium makes its way back to t he center where it can crush the PWN. This late r 



phase of evol ution has recently b een the subject of detailed studies ( van der Sw aluw et al.. 200 1Blond in et al.. 2 001). 
In particular, iBlondin et al. (200lh noted that the reverse shock is likely to be asymmetric so that the PWN can be 
displaced from its position over the pulsar. This scen ario provides an explanation for the displacement of the radio 
emitting PWN in the Vela remnant feock et al.. 1 998) and other remnants. 

Here, I emphasize PWNe that are likely to be in the earlier phase of evolution, before the reverse shock effects. 
Recent discoveries at X-ray and radio wavelengths have substantially increased the number of such objects. In § 2, 
the possible members of this class are listed. For objects with an approximately constant pulsar power, the expansion 
in a supernova is treated in § 3. Constraints implied by the energy in the nebulae are discussed in § 4. The conclusions 
are in § 5. 

YOUNG PULSAR WIND NEBULAE 



A list of probable young PWNe in which central pulsars have been identified is given in Table 1 ; these objects 
are plausibly interacting with ejecta. The second column gives the observed pulsar period, P, and the third column 
gives the characteristic pulsar age, t c h = P/2P. If the pulsar is born rotating much more rapidly than the current 
rate and the braking index is n = 3, then t c h is the actual age. If the pulsar is born with a period close to its current 
period, it can be younger than t^. Alternatively, if the pulsar is born spinning rapidly and has a braki ng index n < 3 , 
it can be older th an t r h. The fourth colu mn is an estimate of the actual age. For 0540-69 dKirshner et al.. 1989) 
and G292.0+1.8 (IMurdin and Clark. 19791) . the age estimate is from the size of the nebula and velocities found from 
optical spectroscopy. For the other objects with ages, the estimate is from a tentative supernova identification, given 



Table 1 . Young pulsars/pulsar wind nebulae 



Object 


P 


P/2P 


Age 


SN 


SNR 


Swept up 


Refs. 




(msec) 


(year) 


(year) 






ejecta? 




0540-69 


50 


1660 


760 




Yes 


Yes (optical) 


1 


3C 58 


66 


5390 


821 


1181 




Yes (X-ray) 


2 


Crab 


33 


1240 


948 


1054 




Yes (optical) 


3 


Kes75 


325 


723 






Yes 




4 


G292.0+1.8 


135 


2890 


^1600 




Yes 


Maybe (optical) 


5,6 


Gll.2-0.3 


65 


24,000 


1616 


386 


Yes 




7 


MSH 15-52 


150 


1700 


1817 


185 


Yes 




8 


G54. 1+0.3 


137 


2890 










9,10 



References: (1) IKirshner et al. (1989h: (2) iBocchino et al. (200lh: (3 ) ISankrit and Hester (1997 ): 
elfand et al. (2 0oi: (5) IC amilo et al. (2002hh: (6) IHughes et al. (2001 h : m iRoberts et al. (2002h : 
Gaensler et al. (200 lh : (9) Lu et al. (2002h : (lO)ICamilo et al. (2002*]) 



(4) 
(8) 



in the next column. The identification of the Crab with SN 1054 is generally considered very secure, but the other 
identifications are less secure. There is stil l some uncertainty over whether all of the events are in fact supernovae, 
e.g., SN 386 (Stephens on and Green. 20o'2T) . The next column indicates whether an extended supernova remnant is 



observed around the PWN, and the penultimate column indicates whether there is evidence for ejecta swept up by the 
PWN. 

There are indications that these nebulae are in the early phase of interaction with freely ex panding ejecta. One ex- 

pectation in this picture is that the PWN should shock and sweep up a shell of supernova ejecta jRevnolds and Chevalier. 1984Cheva 
The shock wave is initially expected to be radiative, but it becomes nonradiative due to the decline in the super- 
nova density. A good example of an ejecta shell is the complex of f ilaments in the Crab Nebula. In this case, 
there is evidence for the shock wave in the ejecta JSankrit and Hester. 1997h : it may be in the process of making 
a transition from radia tive to nonradiative. For 3C 58, there is evidence for X-ray emission from swept-up ejecta 



(Bocchin o et al.. 20011) . imp lying the presenc e of a no nradiati ve shock. Th e optical emission from 0540-69 may be 



from a radiative shock (IKirshner et al.. 198^Chevalier and Fransson. 1992h. The optical e mission from G292.0+1.8 
appears to be from the vicinity of the PWN (M urdin and Clark. 19 79 Hugh es et al.. 20011) . but the relation between 
them has not yet been determined. 

Another expectation of the model with expansion in ejecta is that the pulsar should be centrally located within 
the PWN. This appears to be true for the objects listed in Table 1, although in most cases the PWN has an asymmetric 
boundary. The PWN in G292.0+1.8 is substantially off the center of the surrounding s upernova remnant, which has 



led to the suggestion that the pulsar has a velocity ~ 770 km s _1 (IHughes et al.. 200 lh . If this is the case, the pulsar 
is expected to move to a place in the ejecta where it is comoving with the surrounding gas and it is surrounded by 
uniformly expanding ejecta. Some degree of asymmetry may develop if there is a gradient in the surrounding density 
distribution. 

MODELS FOR 3C 58 AND Gll.2-0.3 

If we assume that 3C 58 and Gll.2-0.3 are associated with SN 1181 and SN 386, respectively, models for the 
PWNe can make use of the fact that the characteristic age is much larger than the true age, so that the pulsars have not 
significantly spun down. This allows the assumption that the pulsar power, E, is constant during the evolution. 

We begin by considering the expansion of the PWN into the inner parts of the sup ernova ejecta. The density dis- 
tribution into which the PWN initially expands can be estimated from explosion models. Chevalier and Fransson (1992h 



used simple power law models for the density distribution. More detailed models have been considered bv lMatzner and McKee (199 
who give asymptotic forms for the the final density distribution at low and high velocities. For the cases to be studied 
here, the asymptotic low velocity density profile is applicable over the time of interest. For an explosion in a star with 



a radiative envelope, the inner density profile, using eq. (46) of 



Matzner and McKee (1999), can be expressed as 



x -1.06 / M . \ 1.97 

pt? = 4.3 x10 s =- E^gcm-^s 3 , (1) 

where v = r/t is the free expansion velocity, M e j is the total ejected mass, and £51 is the explosion energy in units of 
10 51 ergs. For a star with a convective envelope, the density distribution may be flatter, but the density at 1000 km s _1 
is close to the above value. 

The expansion of a PWN in a density distributi on can be approximately treated as a thin shell driven by a uniform 
pressure wind bubble with adiabatic index 7 = 4/3 ( Ostriker and Gunn. 197llChevalier. 197'feevnolds and Chevalier. 1 9841 



The radius can be found analytically if the pulsar power, E, is constant and the surrounding medium has a power law 
density distribution. As discu ssed above, this may be the case for 3C58 and Gl 1.2-0.3. Then, from eq. (2.6) of 
IChevalier and Fransson (19 92). we have 

R,-0X>BZ»1&>(J* S -) (^) PC (2) 

where E^s is E in units of 10 38 erg s _1 . In cases where R p , E, and t are known, we can solve for M e j / '-Efj 45 '; 
because E51 ~ 1 is expected, we have an estimate of the total ejecta mass. These quantities are known for 3C 58 and 
Gl 1.2-0.3 if the historical supernova identifications are assumed, and the results are given in Table 2. These estimates 
involve a spherical approximation for P WNe that are apparently not spherical, but an a verage radius can be taken 
(Wolt ier et al. . 1997 Robert s et al.. 2Q02K The assumed distances are 3.2 kpc for 3C 58 (Robert s et al.. 19931) and 5 



kpc for Gl 1.2-0.3 dGreen et al.. 19881) . The expected values of M ej for a core collapse superno va are typically severa l 
M & or more. In the case of a very fast Type Ic supernova, SN 19941, M e j was only ~ 1-M© Jlwamoto et al.. 1994T) . 
which is an extreme case. The point is that the model leads to a value of M e j that is smaller than expected for 3C 58. 
The radius is larger than would be expected if it were expanding into a normal supernova. There is no problem with 
Gl 1.2-0.3. 



Table 2. Estimated ejecta and swept up mass 



Object 


-^38 


Radio radius 


M ej /E$f 


Predicted 


Observed 






(pc) 


(M ) 


M sw (M ) 


M sw (M ) 


3C 58 


0.27 


3.3 


0.1 


0.002 


0.1 


Gl 1.2-0.3 


0.064 


0.9 


3.5 


0.05 





Another constraint comes from the amount of mass swept up by the wind bubble, M sw . An integration over the 
central density shows that M sw l.0ER~ 2 t 3 , fairly independent of supernova density distribution. Values of M sw 
deduced in this way for 3C 58 and Gl 1.2-0.3 are given in Table 2. The value of 0.002M^ deduced for 3C 58 can 
be compared to the O.1M found from X-ray observations (IBocchino et al.. 20011) . Again, there is a problem with 
the model mass being too low. The mass could be brought into agreement with the observed value if the age were 
increased by a factor ~ 3. 

These models assume that the supernova ejecta are swept into a thin shell that remains spherically symmetric. 
However, the shell is being accelerated and is subjected to Rayleigh-Taylor instabilities, which can decrease the 
coupling between the pulsar wind bubble and the swept up gas. In the limit that there is no further acceleration after 
the ejecta are shocked, the PWN radius (eq. [2]) is increased by a factor of 1.4. The masses M e j and M sw are 
increased by a factor of 2, which does not change the conclusion about the difficulties with 3C 58. 

The model with constant E can be used to predict the internal energy in the PWN. This energy is reduced from the 
total deposited energy, Et, because of work done on the surrounding supernova gas. For a range of flat central density 
distributions, the internal energy is 0A5Et (Table 1 of Chevalier and Fransson 1992). The internal energy in a PWN 
has relativistic particle and magnetic field components; a minimum value for the total energy in p articles and fields 
can be found from the synchrotron luminosity and the emitting volume (e.g., Pacholczyk 1970). Tarn et al. ( 2002) 



discuss the radio emission from the PWN in Gl 1.2-0.3. The value of the minimum energy deduced from the radio 
emission for the two PWNe is given in Table 3. The actual energy must be larger when a larger frequency range is 
considered. It can be seen that there is not a problem with the energy for Gl 1.2-0.3, but that the energy in 3C 58 
appears to be larger than that expected for the observed pulsar and the designated age. A larger age for 3C 58 would 
allow a larger energy to be deposited in the nebula. 



Table 3. Minimum energy in radio emitting particles (10 — 10 Hz) 



Object 


Lradio 


Et 


E ■ 


■^min/ Et 




(10 34 ergs s- 1 ) 


(10 48 ergs) 


(10 48 ergs) 




3C 58 


2.8 


0.7 


1.0 


1.5 


Gl 1.2-0.3 


0.12 


0.3 


0.03 


0.1 



There are thus several arguments for 3C 58 being older than SN 1181, even though Stephe nson and Green (2 002) 



consider the identification to be secure. The problems are that the PWN is too large to be expanding into a normal 
supernova, the expected mass swept up by PWN smaller than observed, and the internal energy is larger than can be 
suppli ed by the pulsar. A la rger age for th e remnant is consis tent with the slow expansion of 3C 58 observed at both 
radio (Bietenh olz et al.. 200 ll) and optical ( Fesen et al.. 1988h wavelengths. 



ENERGY EQUIPARTITION IN PWNe 

One of the important properties of a PWN is the relative amount of energy in particles and in magnetic fields. 
In the Crab Nebula, there is approximate equipartition overall in these energies. In the detailed MHD model of 



Kennel and Coroniti (1984) for the Crab Nebula, this property is produced by the choice of the a parameter, the ratio 
of Poynting flux to particle kinetic energy flux in the pulsar wind. The value a « 0.003 is deduced; the magnetic field 
is relatively weak in the wind and is increased by the shock compression and further compression in the decelerating 
postshock flow. This value of a is close to the upper limit that is allowed in this kind of model, or the flow would 
not be able to decelerate to meet the outer boundary condition, but there is no particular reason for this value to be 
produced. 

One way to estimate the overall magnetic field in a PWN is from the synchrotron break frequency, z^ r , and the 
age of the PWN. The determination of depends on the spectrum of the nebula. There is increasing evidence that 
the particle spectrum injected into PWNe typically has at least one intrinsi c spectral break. Models for th e radio to 



X-ray emission with a single power law injection spectrum generally fail (Reynolds and Chevalier. 1984), and the 



well-observed Crab Nebula spectrum requires an injection spectrum with a break (e.g., Amato et al., 2000). With 
synchrotron losses, the spectrum develops a further break. 

Two PWNe for which there is age information and information on z^ r are 3C 58 and MSH 15-52. As discussed 
above, 3C 58 may be older than 821 years, which I consider as a lower limit; a lower limit on the age yields an upper 
limit on the magnetic field. MSH 15-52 has been sugges ted to be the remnant of SN 185, although G315.4-2.3 is 



another candidate for the remnant of this supernova (Stephenson and Green, 2002). However, the large size of the 
nebula associated with MSH 15-52 (Gaensler et al.. 200 lj) indicates that it is not significantly younger than 1700 



years. 

The determination of vy. depends on the interpretation of the overall spectrum of the PWN. For 3C 58 and MSH 
15-52, as for most PWNe, there are detections at only radio and X-ray wavelengths. However, the fact that the extent 
of the X-ray emission is comparable to that of the radio emission in both cases indicates that iv is not much lower 
than X-ray energies, assuming that particles originate close to the pulsar and move out in the nebula. In the case of 
3C 58, the X-ray spectrum steepens at large radii, showing that synchrotron losses are significant in the X-ray regime 
dTorii et al.. 20 00). An estimate of hvbr is thus ~ 0.5 keV, which is consistent with the fact that the X-ray spectrum 
is somewhat steeper than the spectrum from radio to X-ray wavelengths. The implications for the magnetic field 
and magne tic energy are given i n Table 4. The magnetic field strength in MSH 15-52 can be estimated from similar 



arguments dGaensler et al.. 2001 ) and is given in Table 4. In this scenario, the low frequency breaks are intrinsic to the 



injected particle spectrum; this is a controversial point and there have been discussions of the low frequency breaks in 



terms of synchrotron losses (e.g., Woltjer et al., 1997 on 3C 58 and Roberts et al., 2002 on Gl 1.2-0.3). 

Table 4. Comparison of magnetic and minimum internal energies 



Object hv hT B E B E min 
(keV) (fiG) (10 47 ergs) (10 47 ergs) 

3C58 05 16 04 10 

MSH 15-52 1 8 0.3 7 



It is also possible to estimate minimum magnetic plus particle energy, E m i n , from the radio synchrotron emission 
(10 7 -10 n Hz). The results, given in Table 4, show that the magnetic energy is considera bly less than the total internal 
energy in both cases, so that the nebulae are particle dominated. In the context of the Kennel and Coroniti (1984) 



model, this would requ ire a remarkably low Poynt ing energy flux in the pulsar wind. However, the result is consistent 



with the finding in the Kennel and Coroniti (1984) model that a particle dominated wind is needed to produce a shock 



front and deceleration of the flow to match the outer boundary. This suggests that pulsar winds may have a range of 
magnetizations, giving rise to a range of nebular properties. If cases with a highly magnetized wind occur, they would 
give rise t o something other than a standard pulsar wind nebula. In this limit, the wind termination shock moves in to 
the pulsar ( Rees and Gunn. 1974Emmering and Chevalier. 1987) so that the immediate surroundings of the pulsar are 
in communication with the ambient medium. 

DISCUSSION AND CONCLUSIONS 

The increasing number of PWNe are providing many examples which can be compared to models initially devel- 
oped for the Crab Nebula. In a model with interaction with freely expanding supernova ejecta, the PWN properties can 
provide a check on the age esti mate for the nebula. Such mod els have previously been developed for the Crab Neb- 
ula, MSH 15-52, and 0540-69 Jchevalier and Fransson. 1992h . Recent observational results on 3C 58 and Gl 1.2-0.3 
allow similar models to be considered for these objects, which are of special interest because of their possible iden- 
tifications with historical supernovae. The models indicate that 3C 58 is older than SN 1181, but that G 11.2-0.3 is 
consistent with being the remnant of SN 386. These tentative conclusions need to be followed up by more detailed 
studies of the remnants. In both cases, the model predicts that the P WN is driving a sh ock front into freely expanding 
ejecta. Gas shocked in this way may have been observed in 3C 58 dBocchino et al.. 200 lh . but further observations 
are needed. Such gas has not yet been observed in Gll.2-0.3. The external s upernova remnant in teraction also pro- 



vides constraints on the system. This interaction is clearly seen in Gll.2-0.3 (Rober ts et al.. 2 002). but not in 3C 58, 
implying interaction with low density surroundings for that case. 

Another finding here is that the 3C 58 and MSH 15-52 PWNe are particle dominated, which requires that the 
pulsar winds have a very low magnetization parameter. This property may be one of t he reasons w hy these PWNe have 



low efficiencies of X-ray luminosity production compared to the pulsar power ( Chevalier. 2000h . 3C 58, G292.0+1.8, 



Gll.2-0.3, MSH 15-52, and G54.1+0.3 all have significantly lower levels of X-ray luminosity production efficiency 
than the Crab Nebula. 
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